Fano resonances and Aharonov-Bohm effects in transport through 

a square quantum dot molecule 

Z. Y. Zeng 1 ' 2 , F. Claro 1 , and Alejandro Perez 1 
l.Facultad de Fisica, Pontificia Universidad Catolica de Chile, Casilla 306, Santiago 22, Chile 
2. Department of Physics, Hunan Normal University, Changsha ^10081, China and 
CC AST (World Laboratory). P. O. Box 8730, Beijing 100080, China 

(February 1, 2008) 

Abstract 

We study the Aharonov-Bohm effect in a coupled 2x2 quantum dot array 
with two-terminals. A striking conductance dip arising from the Fano inter- 
ference is found as the energy levels of the intermediate dots are mismatched, 
which is lifted in the presence of a magnetic flux. A novel five peak structure 
is observed in the conductance for large mismatch. The Aharonov-Bohm evo- 
lution of the linear conductance strongly depends on the configuration of dot 
levels and interdot and dot-lead coupling strengths. In addition, the magnetic 
flux and asymmetry between dot-lead couplings can induce the splitting and 

combination of the conductance peak(s). 
PACS numbers: 73.23.-b,73.63.-b,73.63.Kv 
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I. INTRODUCTION 



The famous effect predicted by Aharonov and Bohm (AB) in 1959 Q concerns the es- 
sential influence of a vector potential on the interference pattern of two beams of electrons 
confined in multiply connected paths, over which the magnetic field is zero. It is manifested 
by the periodic oscillation of the conductance of a ring or a cylinder as a function of the 
enclosed magnetic flux $ . It has been shown that the AB effect also exists in some singly 
connected geometry such as a point contact or a disk shape in a two-dimensional electron 
gas, in which circulating edge states enclose a well-defined magnetic flux [^J. A significant 
hallmark in mesoscopic experiments is the phase measurements of the transmission ampli- 
tude through a quantum dot embedded in an AB ring in the Coulomb blockade regime [|J . 
It revived interest on the AB effect in condensed matter physics ||. 

Quantum dots are highly tunable artificial mesoscopic structures ||, usually called ar- 
tificial atoms or artificial molecules. Some novel physical properties, such as the Coulomb 
blockade and the turnstile effect, have been demonstrated in some elegant experiments ||. 
At low temperatures, transport through a single quantum dot is dominated by resonant tun- 
neling and Coulomb blockade. When the condition for resonant tunneling is not accessible, 
co-tunneling events involving the simultaneous tunneling of two or more electrons through a 
virtual intermediate level can be dominant 0. When quantum dots are connected in series 
with tunneling barrier (s), it is expected that the coupling between dots plays a crucial role 
in the transport properties of the coupled structure. On the other hand, electron tunneling 
through the coupled dot system is very sensitive to the incoherent scattering, which has 
trivial effects on the transport properties of a single dot [§]. In the past decade, transport 
through coupled quantum dot systems has received considerable amount of attention |9|-|lT| . 

In this paper, we investigate coherent electronic transport through four coupled quantum 
dots located at the four corners of a square enclosing a magnetic flux. This kind of coupled 
dot structure is the basic unit of two-dimensional quantum dot arrays. Our goal is to find 
out how the conductance of such a structure depends on the coupling strength between dots, 
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the arrangement of dot levels and the magnetic flux. In the present work, we do not take 
into consideration the interdot and intradot electron-electron interactions. We also assume 
that just one level is relevant on each dot. This allows us, on the one hand, to obtain some 
exact analytical results. On the other hand, we can gain much more clear physical insight 
into the dependence of the conductance of the system on the structure parameters (interdot 
couplings and dot levels) and the magnetic flux. Though the model we consider is simple, 
some novel and interesting results are obtained arising from the interplay of the specific 
configuration of dot levels and interdot couplings, and the magnetic flux. 

The rest of the paper is organized as follows. In Sec. II we derive the current through the 
2x2 quantum dot array using the Keldysh nonequilibrium Green's function formalism and a 
recursive Green's function technique. Section III presents the linear conductance spectra of 
the dot array in the symmetric and asymmetric coupling case. We also investigate in detail 
the Ahoronov-Bohm oscillation of the linear conductance in some specific configurations of 
the dot levels and the interdot couplings. Our concluding remarks are given in Sec. IV. 



II. FORMULATION 



We consider four quantum dots in a square array |12| enclosing a magnetic flux $ (Fig. 
1), with two dots facing each other connected to the left and right leads. For simplicity, we 
ignore the intra- and inter-dot Coulomb interactions, and assume that just one energy level 
is relevant at each dot. Then the hamiltonian describing such a system is 



H — H Dot + H Lead + H Dot 

—Leadi 

4 

H Do t = e idfdi + (ti 2 dfd 2 + t^e^di d 3 + t 24 djd 2 + £346^3 + h.c), 



i=i 

HLead = X! € k a k a k + X! e pK ^Pi 
k&L p€R 

Hoot —Lead 

fcgL peR 
where is the energy level of dot i, ak(a k ) ,b p (bp) and di(df) {i = 1,2,3,4) are the annihi- 
lation (creation) operators of electrons in the left(L) lead, right(R) lead and quantum dot i, 
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and Uj and Vk\ (W P 4) are the inter-dot and dot-lead coupling matrix elements, respectively. 
In Eq. (1), a factored (0 = 27r$/$ 0) = h/e) is attached to £i 3 to account for the magnetic 
flux $ through the structure. 



The current flowing from the left lead to the structure is |L3 
Hp r df 

J L = -f-J ^L(e){f L (e)[G r dldi (e) - G a didi (e)] + G< di (e)}. (2) 

The notation we use throughout is G r £ a Y (t,t') = ^i9(±t =p t') < {X(t), Y + (t')} >, 
G< Y (t,t') = % < Y+(t')X(t) >,X,Y = a k} b P ,d l (z = 1,2,3,4), T L (e) = 27r\V kl \ 2 5(e - e fc ), 
while = [1 + e^ _Mi ^ //fcsT ] _1 is the Fermi-Dirac distribution function of the left lead, 

with chemical potential /i^. The factor of 2 in Eq. (2) is due to the spin degeneracy. 

Following the recursive decoupling technique developed in Ref. [fL4H , the retarded Green's 
function appearing in Eq. (2) G r dldl (e) can be calculated from Dyson's equation 

G r dldl = [e-e l -Z r dl }-\ (3) 
54 = ~T L + \t 12 \ 2 G d2d2 + \t 13 \ 2 G d3d3 + t 12 tl 3 e^G r d2d3 + ty 13 e^G r d3dr (4) 

Here and in what follows G 1 ^. denote the Green's function which is decoupled from the 
quantum dot k (k = min(i,j) — 1). In the above and following equations we drop the 
argument e, and recover it wherever necessary. Following a similar recursive decoupling 
method, one has 

^ r = (e-e 3 )(e-e 4 + |r fl )-t34t* 4 

**> ( e _ e2 )( e _ e3 )( e _ e4 + i ri? )-t 24 t* 4 ( e - e 3)-t3 4 t* 4 ( e - e2 )' 1 > 

q T = (e - e 2 )(e - e 4 + |Fr) - , . 

d3d3 (e - e 2 )(e - e 3 )(e - e 4 + f T R ) - t 2A t* 24 {e - e 3 ) - t u t* u (e - e 2 ) ' U 

** ( e _ e2 )( e _ e3 )( e _ e4 + ir R )-t 24 t* 4 ( e - e3 )-t3 4 t* 4 ( e - e2 )' 1 ) 

qt _ ^34^24 /g\ 

( e _ e2 )( e _ e3 )( e _ e4 + i ril )-t 24 t* 4 ( e - e 3)-t3 4 t| 4 ( e - e2 )' 1 ' 

where T^(e) = 27r|W / P 4| 2 5(e — e p ). In the following we will assume that the coupling matrix 
elements (V^i, W P 4 and ty) are real since the tunnel rates just depend on the amplitude of 
the coupling matrix elements. Substitution of Eqs. (5)-(8) into Eq. (4) yields 



yr , v = _ i_ r (e - e 4 + |r R )[^ 2 (e - e 3 ) + t 2 n (e - e 2 )] - (t 2 2 t 2 u + t 2 3 t 2 4 ) + 2cos(j)t 12 t 13 t 24 t 34 

dl[e) ~ 2 L+ ( e _ e2 )( e _ e3 )( e _ e4 + ir R )-4 (e _ e3) _ t 2 4(e _ e2) 

(9) 

It can be seen from the above expression that the AB effect manifests itself in an additional 
term 2(coscj) — 1)^12^13^24^34 in the self-energy of dot 1 Green's function, the dominant 
contribution to the spectral density of the structure. 

We now have to calculate the lesser Green's function G d dl . It can be obtained from the 
Keldysh formula 

G< (~tr y< /~ia sp< \/~ir |2 f~\r\\ 

didi ~ ( - T didi Zj di LT didi — ^di I ^didi I > \ LV ) 

where 

E< = if L (e)T L + £ tlG< di + t X2 t x ^G< 2d3 + t 12 t 13 e^G< d2 . (11) 

i=2,3 

The lesser Green's functions G d . d] = 2, 3) can be calculated from the following Keldysh 
equation 

Gdidj = X! G didm Y, dmdn G r dnd] , i,j = 2,3 (12) 

m,n=2,3 



where 



f%* = ; fR ^l Rt ^ . m ,n = 2,3 (13) 
(< - < ,) 2 • I 7,, 1 



Here /#(e) = [1 + e( e ~ m ^' kBT ]~ 1 is the Fermi-Dirac distribution function of the right lead, 
with chemical potential //r. It is worth noticing that the above procedure of calculating the 
various kinds of Green's functions can be verified to give the same results as the equation- 



of-motion method [15|. Combining Equations (6) — (10) we find 



GU = (ih(e)T L - 2if R (e)Im[E di + ^T L ])\CT dldl | 2 . (14) 



Notice that 

1 



- no. nr 

^dxdx ^did! 



G r dldl - G a dldl = - j^-W^ = 2Um[^ di }\G dld f. (15) 
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Substituting Eqs. (14) and (15) into Eq. (2), one then obtains the current entering the 
structure from the left lead , 

J=J t = 2 ij d <f^) - fRteWte), (16) 

in which 

r L (e)r R (e) 

T 6 = F ( \ T F ? ^ 6 ' 

r L (e) + r B (e) 

2(r L (e) + r B (e)) r , , z , .,, 1 |2 

is the transmission probability for an electron passing through the quantum dot square 
structure. Equations (16) and (17) are the central results of this work. It can be seen 
that the transmission probability and current are proportional to the dot-lead couplings 
(r^lV) and a generalized spectral density function p(e) of the system. Note that the 
current includes contributions from the upper arm path (1 — > 2 — > 4), and the bottom arm 
path ((1 — ► 3 — > 4), which interfere with each other. At zero temperature, the conductance 
in the linear regime limit ^ ~ /iR — > e F becomes 

h T L + T R 



In obtaining Eq.(18) the wide-bandwidth approximation |13| has been used, i.e., the 
linewidths and Tr are constants independent of energy. From the expression of the 
self-energy (9) one can find that the linear conductance {?($) is a periodic function of the 
magnetic flux with periodicity $o — h/e, and = G {—*&)■ These observations are 

consistence with the general results of the two-terminal setup of the AB ring ||. 

When quantum dot 2 ( or 3) is decoupled from the system, i.e., ti 2 = t 2 4 = ( or 
tn = ^34 = 0), we recover the results for the coupled triple quantum dot chain [|14[ 

J = ^ = x / *Lfr(«o - (19) 

where the transmission probability is 

%(e) = T L T R \t u \% 4 \ 2 A(e). i = 2 {or 3) (20) 



Here 



M ^) 2 + r|/4]/{[(6 - eQ[(e - e 4 ) 2 + I*/4] - t 2 4 (e - e 4 )] 2 + |t^| 4 r^/4} 

1 J " (e - e x - t 2 u ReB(e) 2 + (T L - 2t 2 l ImB(e) 2 / '4 

i = 2 (or 3) (21) 

where 

B{e) = {e-e i ^-j' 1 ■ * = 2 (or 3) (22) 

When the levels of the three dots are aligned, we find the condition for complete transmission 
(T 3 = 1) to be 

| = . z = 2 (or 3) (23) 

It can be expected that the criteria for complete resonant transmission or maximum peak 
conductance plays an important role in the formation of conductance peaks with the con- 
ductance quanta 2e 2 /h for the square quantum dot array. 

III. LINEAR CONDUCTANCE CHARACTERISTICS 

A. Symmetric coupling between dots and leads 

It is well known that the two-terminal conductance of an AB ring is a periodic even 
function of the enclosed magnetic flux $ p. The oscillation of the conductance, with 
periodicity h/e, results from the quantum interference between the paths located in the upper 
and lower arms of the ring. For electron tunneling through a single dot, the conductance 
peak reaches its maximum when = Tr. In the coupled double quantum dot structure 
with aligned energy levels, the conductance peak reach a maximum value and splits into a 
double peak with peak separation 2{t 2 — T 2 /4) 1//2 when the interdot coupling t becomes larger 
than the width of the peak T in the symmetric dot-lead coupling case [0. For the three 



coupled dot string, the situation is more complicated, though the condition for maximum 
conductance peak is the simple Eq. (23) since the levels of the dots are aligned. We expect 



that the conductance peak will split into three maxima when the interdot coupling becomes 
larger than the width of the peak, as in the double dot case. 

For further comparison, we first calculate the linear conductance of the quantum dot 
array in the absence of the magnetic flux $ and in the symmetric dot-lead coupling case 
(T L = T R ). Throughout this paper, energy and coupling strengths are measured in units of 

= T. The results for aligned levels and various kinds of interdot coupling configurations 
are shown in Fig. 2 (a), while Fig. 2 (b,c) are for different arrangements of the dot levels. 
When the four dot levels are aligned, three conductance peaks appear if ti 2 = ti3 and 
^24 = ^34- Then the conductance possesses the same three resonant peak structure as the 
three dot string. One can also find that the criteria for complete resonant transmission for 
the three-dot systems manifests itself in our four dot structure. This is not strange since the 
conductance in our case is simply the superposition of the contributions from the upper and 
lower three-dot arms. When tyi = £34 < ^13 = £24, the situation becomes different. In this 
case, one can consider the dot structure as a system of two coupled quantum dot dimers. 
One observes in the conductance spectra four peaks, arranged in two groups separated by 
a conductance gap. The separation between the two group centers is mainly determined by 
the larger interdot coupling, and the peak separation of each group by the smaller one. The 
peak and group center separations have trivial dependence on the dot-lead couplings Tl and 
T R in the coupled quantum dimers case. The conductance gap is introduced by the larger 
coupling strength between the two dimers and can not be lifted by a magnetic flux, contrary 
to the case of the Fano conductance dip discussed below. The gap can be broadened by 
increasing the coupling strength(s) ti 2 ( and/or t 2 i) between the two dimers. 

Figures 2 (b,c), displaying cases of mismatched levels, show more interesting conductance 
spectra. When t 12 = £13 = t 2 4 = Ha a striking novel conductance dip {Q = 0) appears 
halfway between the energy levels of dots 2 and 3 (Fig. 2(c) and dashed line in 2(b)) if 
these levels are different. When tyi = t 2 4 « ti 3 = t 34 the conductance dip is pinned 
at the energy level of the quantum dot with small couplings to the system (solid line in 
Fig. 2 (b)). The middle peak appears split into two asymmetric peaks. The asymmetry 
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between the two split peaks depends on the mismatch between the energy levels of dots 
2 and 3, as well as the arrangement of the levels of other dots 1 and 4. Novel five spikes 
( the dashed line in Fig. 2 (b)) in the conductance can be seen if the mismatch between 



the levels of dot 2 and 3 is large enough. The dip results from a Fano interference [16 



between two distinct current paths, one being direct nonresonant and the other nondirect 



resonant fll7|Jl8|| . The conductance can be generally written fl7j as Q = Qnon ^-tl] , where 



Qnon is the nonresonant conductance, e = e — e reso , and q is the 'Fano parameter' assessing 
the asymmetry of the lineshape. Two limiting cases exist [Ref. 17(b)]: (a) q — > oo, with 
a dominating resonant transmission resulting in Breit-Wigner resonances and, (b) q — > 0, 
with a dominating nonresonant transmission leading to asymmetric dips. Since here t\2 = 
tn = ^24 = £34, one can consider either the arm consisting of dot 2 or 3 as the resonant 
(nonresonant) path. This special symmetry locates the Fano-type conductance dip at neither 
the level of dot 2 nor that of dot 3, but halfway between them. On the other hand, when 
t\2 = tiz = ^24 = ^34, we find from Eqs. (9), (17) and (18) that Q oc ^l^"^"! 3 ^) • It is evident 
that two resonance peaks exist when ep = £2 and €f = £3, while Q=0 when €f = (62 + e 3)/2, 
which is just the position of the Fano conductance dip. For large mismatch between energy 
levels of dots 2 and 3, Fano interference between the two paths develops both a conductance 
dip and a peak halfway between these levels 62 and £3. One can then observe a five peak 
structure in the conductance. When the levels of dots 2 and 3 are matched, the resonant 
peak structure of the dot systems behaves to a large extent the same way as the coupled 
three-dot system [ 14] . 



Now we turn to study the effect of the magnetic flux $ on the conductance spectra of 
our square dot system. Figure 3 exhibits clearly the AB oscillations of the conductance 
as a function of the magnetic flux <fi within a period at ep = 0.1T. When the two arms 
of the dot ring are symmetric, i.e., the dot levels and interdot couplings are arranged to 
be the same, the conductance Q vanishes at <p = ir (Fig. 3(a) and the solid line in Fig. 
3(c)). In the case of asymmetric arms, the conductance is nonzero within the whole period. 
The vanishing of Q at = it results from the complete destructive interference between 
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the two symmetric arms, which will be analyzed in more detail below. If the two arms are 
not symmetric, complete destructive interference is never achievable, and the conductance 
is always nonzero. The more asymmetry between the two arms, the larger the conductance 
near = ir. To get a clearer insight into the influence of a magnetic flux on Q, we give 
in Figs. 4 and 5 the conductance of the structure as a function of the Fermi energy for 
different values of magnetic flux within half a period (0 = 0, 7r/4, 37r/4, it), and different 
configurations of interdot couplings and dot levels. The dependence of the conductance on 
the magnetic flux <3> is apparent. The magnitude of the AB oscillation strongly depends on 
the matching condition of the dot levels and the configuration of the interdot couplings. It 
also differs for different pinning positions of the Fermi energy at the leads. The presence 
of a magnetic flux can lead to the splitting (Fig. 5(b)) and combination (Fig. 4(c) and 
Figs. 5 (a,c,d)) of conductance peak(s) since it changes the interference between the two 
current paths, as determined by term 2(cos0 — 1)^12^13^24^34 in Eq. (9). Another striking 
effect induced by the magnetic flux is that the central conductance peak, if present (Fig. 
4), becomes a conductance dip when the levels of dots 2 and 3 are aligned, while the Fano- 
type dip appearing in the case of mismatched levels is lifted (Fig. 5), contrary to the 
usual conductance gap which survives in the presence of a magnetic flux (Fig. 4 (c)). This 
feature is typically a result of the modified Fano phase interference by the magnetic flux. 
When = 7r, the conductance is depressed everywhere if e 2 = e 3 and £12^24 = ^13^34, as 
shown in Fig. 4 and Fig. 5 (b). The overall depression of the conductance results from 
the complete destructive interference between the upper-arm and lower-arm current paths. 
When e 2 = e 3 = e and £12^24 = ^i3^34 5 one finds Q oc 2(1 + cos^t^t^, which implies zero 
conductance for any value of the Fermi energy if = re. We can also analyze it from the 
interference pattern of the two current paths containing either dot 2 or 3. The respective 
transmission amplitudes T 2 = t^t^/i^F — e 2 + ^0 + ) and T 3 = ti 3 t u e 1 ^ / (ep — e 3 + i0 + ) 
equal Tj{ep — e + i0 + ) and Te 1 ^ /(e^ — e + i0 + ), respectively, when £12^24 = ^13^34 = T and 
£2 — e 3 — e. Then the total transmission probability T = |T 2 + T 3 | 2 for electron tunneling 
through the dot structure is proportional to 2T 2 (1 + cos0)5(ei? — e), which is zero when 
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(j) = 7T. 



B. Asymmetric coupling between dots and leads 

In the preceding subsection we have reported in detail the linear conductance spectra in 
the absence and presence of a magnetic flux, in the limit of symmetric dot-lead couplings. 
There we saw symmetrically- located conductance peaks if the dot levels are aligned, which 
is a consequence of the particle-hole symmetry. In this subsection, we want to study how 
the conductance is modified if the dot-lead couplings are asymmetric. 

We model the asymmetry between the dot-lead coupling I\ and Tr by setting = 
Figures 6 and 7 present the comparison of conductance spectra in the cases of symmetric 
(thin lines) and asymmetric (thick lines) dot-lead couplings, with structure parameters the 
same as in Fig. 4 and 5, respectively. A surprising feature in these figures is that the 
asymmetry between dot-lead couplings can remove the symmetry of the conductance spectra 
around some energy when the energy levels of dots are not aligned (Figs. 7 (c,d)), while 
it can not if they are aligned (Fig. 6). Also depending on the value of the magnetic flux 
0, the conductance in the case of asymmetric dot-lead couplings is lifted (Figs. 6 (a,c), 
4> = 37r/4) or suppressed (Figs. 6 (a,c), = ir/4 and Fig. 7 (a) = 7r/4, 3tc/A), as compared 
to the symmetric case. We also observe another interesting phenomenon. The conductance 
changes in a nonuniform way when the dot-lead couplings become asymmetric. This is 
reflected in the observation that some conductance peaks are lifted while the others are 
suppressed (Fig. 6 (b) and Figs. 7 (c,d)). In addition, we find that the asymmetric dot-lead 
couplings would introduce a splitting of the conductance peak (Fig. 6 (b), Figs. 7 (c,d)). 
It is worth mentioning that the position of the conductance peaks remain nearly unshifted 
when the dot-lead couplings become asymmetric. 

These features can be understood for linear conductance from the analysis of the gener- 
alized spectral function p(e) in Equation (17), which depends on the dot-lead couplings Tl 
and T/j. In the case of a single quantum dot, a difference in the value of the dot-lead cou- 
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plings and T R does not change the position of the resonant transmission or conductance 
peak, but does alter its height. In our case, the dependence on and T R of the spectral 
function or conductance is much more complicated. However, it appears that a difference 
between T L and T R does also modify the weight of the peaks in the generalized spectral 
function p(e) in a complicated way dependent on the applied magnetic flux. The position 
of the conductance peaks is kept intact by the difference between Tl and T R . 

IV. CONCLUSIONS 

To sum up, we have investigated the transport properties in the linear regime, of a 
two-terminal quantum dot square enclosing a magnetic flux. Although our results are not 
directly applicable to more complex arrays involving more conducting paths and dots , we 
believe that our study will provide insight into the physics of interference that such arrays 
will exhibit. Our main results may be summarized as follows: (1) The AB oscillations of the 
conductance strongly depends on the configurations of dot levels and interdot and dot-lead 
couplings, along with the position of the Fermi energy of the leads; (2) a striking conductance 
dip is developed due to the Fano interference when the levels of the intermediate dots are 
mismatched, and a novel five peak structure is observed for large mismatch; (3) the magnetic 
flux and an asymmetric dot-lead coupling can induce the splitting and combination of the 
conductance peak(s); (4) the interplay of a magnetic flux and asymmetry between dot-lead 
coupling can lead to interesting AB oscillations of the linear conductance. If the electron- 
electron interactions are taken into consideration, further structure in the conductance may 
be expected due to the presence of an additional energy scale. Results in the regime in 
which such energy scale is relevant will be reported elsewhere. 
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FIGURES 

FIG. 1. Schematic plot of a 2x2 quantum dot structure enclosing a magnetic flux $ with two 
terminals. 

FIG. 2. Conductance Q as a function of Fermi energy ep in the case of symmetric dot-lead 
couplings (Fl = T^) and in the absence of magnetic flux when the levels of the four dots are 
aligned (a) and not aligned (b, c). 

FIG. 3. AB oscillations of the conductance Q as a function of magnetic flux cp in the case of 
symmetric dot-lead couplings (r^ = Fr). The structure parameters are the same as in Fig. 2 and 
the Fermi energy €f is set to be O.ir. 

FIG. 4. Evolution of the conductance Q with the Fermi energy €f for the case of symmetric 
dot-lead couplings (r^ = Fr) and the levels of the four dots aligned, and = (solid), <p = 7r/4 
(dashed), cp = 3ir/4 (dotted) and (ft = it (dash-dotted). 

FIG. 5. Same as Fig. 4, but with the levels of the four dots not aligned. 

FIG. 6. Conductance Q in the cases of symmetric dot-lead couplings (r^ = Tr, thin line) and 
asymmetric dot-lead couplings (Tr = Tl/2, thick line) for aligned four dot levels. The solid and 
dotted lines correspond to = 7r/4 and 37r/4, respectively. 

FIG. 7. Same as Fig. 6 but for non-aligned levels. 
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Fig. 1 Zeng et al. 



Fig. 3 Zeng et al. 
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Fig. 4 Zeng et al. 
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Fig. 5 Zeng et al. 
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Fig. 6 Zeng et al. 
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Fig. 7 Zeng et al. 




